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Abstract 



In sexually reproducing animals, primordial germ cells (PGCs) are often set aside early in embryogenesis, a strategy that 
minimizes the risk of genomic damage associated with replication and mitosis during the cell cycle. Here, we have used 
germ line markers {piwi, vasa, and nanos) and microinjected cell lineage tracers to show that PGC specification in the 
leech genus Helobdella follows a different scenario: in this hermaphrodite, the male and female PGCs segregate from 
somatic lineages only after more than 20 rounds of zygotic mitosis; the male and female PGCs share the same (meso- 
dermal) cell lineage for 19 rounds of zygotic mitosis. Moreover, while all three markers are expressed in both male and 
female reproductive tissues of the adult, they are expressed differentially between the male and female PGCs of the 
developing embryo: piwi and vasa are expressed preferentially in female PGCs at a time when nanos is expressed 
preferentially in male PGCs. A priori, the delayed segregation of male and female PGCs from somatic tissues and 
from one another increases the probability of mutations affecting both male and female PGCs of a given individual. 
We speculate that this suite of features, combined with a capacity for self-fertilization, may contribute to the dramatically 
rearranged genome of Helobdella robusta relative to other animals. 



Introduction 

The near ubiquity of sexual reproduction among modern 
animal taxa and the involvement of conserved sets of genes 
in specifying their primordial germ cells (PGCs) indicate that 
sexual reproduction via specialized germ cells is ancestral to 
Bilateria. Conserved genes involved in PGC specification in- 
clude homologs of \/asa, piwi, and nanos. Vasa is an ATP- 
dependent RNA DEAD box helicase that plays an essential 
role in the establishment and maintenance of the germ line 
by promoting the translation of germ line-specific mRNAs 
(Carrera et al. 2000; Fujimura and Takamura 2000; Komiya 
et al. 1994). Piwi, characterized by an RNAase H domain, is 
involved in cell autonomous, self-renewing, germ stem cell 
maintenance (Cerutti et al. 2000; Cox et al. 2000; Song et al. 
2004; Seto et al. 2007). Nanos is a translational repressor char- 
acterized by a C-terminal domain comprising two CCHC zinc- 
finger motifs (Wang and Lin 2004; Ewen-Campen et al. 2010). 

Previous work has identified two distinct modes of germ 
cell specification, with names rooted in old debates concern- 
ing the nature of development itself preformation entails the 
segregation of germ line determinants during early embryo- 
genesis while epigenesis entails induction of the germ line at 
later stages by signals from neighboring cells (Extavour and 
Akam 2003). Comparative studies suggest that preformation 
and epigenesis are not mutually exclusive. For example, the 
proposed model of germ line specification for the polychaete 



annelid Platynereis dumerilii entails a combination of in- 
herited cytoplasmic determinants, which establish a mesoder- 
mal posterior growth zone (AAPGZ), followed by inductive 
processes to determine which cells of the AAPGZ become 
the PGCs (Rebscher et al. 2007). Related findings in diverse 
taxa have led to the hypothesis that a broadly conserved 
germline multipotency program (GAAP; involving Piwi, 
Nanos, Vasa, and other gene products) distinguishes a class 
of primordial stem cells (PriSCs), which are intermediate 
between the zygote and the PGCs, and which also give rise 
to somatic cells (Juliano et al. 2010; Solana 2013). 

Leeches of the genus Helobdella (Annelida: Clitellata: 
Hirudinea: Glossiphoniidae) provide another annelid model 
that is evolutionarily distant from Platynereis (Struck et al. 
2011). Clitellate embryos undergo a modified version of un- 
equal spiral cleavage; individual blastomeres can be charac- 
terized reproducibly on the basis of size, position, the pattern 
of cell divisions by which they arise, and their subsequent fates 
(Weisblat and Huang 2001). Intriguingly, genomic analyses 
show extensive genomic rearrangements in Helobdella 
robusta relative to two other lophotrochozoan models 
(Cho et al. 2010; Simakov et al. 2013). 

To begin elucidating mechanisms by which the germ line is 
specified in Helobdella embryos, we had previously identified 
and characterized the expression and function of a nanos 
homolog in H. robusta, finding that after broad early 
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expression, it becomes restricted first to male presumptive 
PGCs and then later appears in female PGCs (Pilon and 
Weisblat 1997; Kang et al. 2002; Agee et al. 2006). Here, we 
have characterized piwi and \/asa homologs in Helobdella: the 
H. robusta genome revealed two copies of each of these genes; 
orthologs were cloned and their expression characterized in 
H. austinensis (whose genes are designated with the prefix 
HoU'); as with Hm-nos, all copies of Hau-piwi and Hau-msa 
are broadly expressed in early stages and gradually became 
restricted to PGCs as development unfolds, consistent with 
the GMP'PriSC model. 

In contrast to our previous characterization o^nams, how- 
ever, expression of p/u// and vasa is first evident in the female 
germ line. Lineage tracing combined with in situ hybridization 
(ISH) revealed that as for the male germ line, female PGCs 
arise from segmental mesoderm. Thus, male and female 
germline fates separate only after 19 rounds of zygotic mito- 
ses, much later than in several other invertebrate models. 
Combined with the capacity of some Helobdella species for 
self-fertilization, our results suggest a mechanism that would 
contribute to the extensive genomic rearrangements seen in 
H. robusta and to the rapid speciation of the genus Helobdella 
relative to other leeches (Oceguera-Figueroa et al. 2011). 

Results 

Sequence Retrieval and Phylogenetic Analyses 
Three pairs of genes encoding homologs of Piwi, Vasa, and 
Argonaut were retrieved from the whole genome sequence of 
H. robusta, and four genes encoding PL10 homologs; we des- 
ignate these genes Hro-piwil, Hro-piwil, Hro-vasal, Hro-\jasa2, 
Hro-agol, Hro-agol, and Hro-plWa, Hro-pllOb, Hro-plWc, and 
HrO'plWd, respectively. Each set of these paralogs was distrib- 
uted across different scaffolds with different flanking genes 
(supplementary table SI, Supplementary Material online); 
this suggests that the gene duplications were not of recent 
origin. To ask whether duplications of the piwi, vasa, and 
argonaut genes had occurred in an ancestor of annelids and 
mollusks, we also explored two other available lophotro- 
chozoan genomes, the mollusk Lottia gigantea and the poly- 
chaete Capitella teleta. Similar to Helobdella, Lottia and 
Capitella, each have two piwi genes. In contrast to Helobdella, 
however, these other two lophotrochozoans each contain 
only one \/asa, one Argonaut, and one plW gene, suggesting 
that these three genes may have undergone duplications in 
the sublineage leading from a last common annelid ancestor 
to leech. 

Consistent with candidate identities, a maximum likeli- 
hood phylogenetic tree groups the piwi and argonaut genes 
in separate clades (fig. 1). The evolutionary dynamics of genes 
within these two clades differ, however. Although the 
arrangement of genes within the argonaut clade suggests 
that multiple, taxon-independent duplications of argonaut 
genes have occurred, the two clear clades separating meta- 
zoan piwil from piwil genes suggests a single, ancestral du- 
plication of the piwi gene (fig. 1A). Similarly, the maximum 
likelihood phylogenetic tree encompassing the \/asa and plW 
homologs groups these two sets of genes into separate clades 



as well. As in the case of the argonaut genes, several duplica- 
tions appear to have occurred independently within both 
these gene clades (fig. IB). 

To study the expression of these piwi- and i/asa-related 
leech genes, we cloned their orthologs from the sibling species 
H. austinensis, which is more easily maintained in laboratory 
culture. Coding sequences of the Hau-piwi and Hau-vasa 
genes were 97-98% similar at the nucleotide level when com- 
pared with their Hro-piwi and Hro-vasa orthologs, respec- 
tively, and 97-98% identical at the amino acid level. 

Expression of Hau-piwil and Hau-vasal in the 
Adult Leech 

Leeches, including Helobdella, are simultaneous hermaphro- 
dites, characteristic of clitellate annelids as a group. Leeches 
have one pair of ovaries (fig. 2A and B), associated with mid- 
body segment 6 (designated as segment M6). The male gon- 
opore lies in segment M5, but the number and distribution of 
the testis sacs varies among (and to a lesser extent within) 
species. Helobdella robusta and H. austinensis normally have 
six pairs of testis sacs, in segments M8-M13 (fig. 2A and B). To 
test the assumption that piwi and vasa are germ cell markers 
in the leech, we processed cryosectioned adult leeches by ISH 
for HaU'Vasal and Hau-piwil. Both genes were expressed 
specifically in the bilaterally paired ovaries and in the testis 
sacs (fig. 2C-F); sense probes gave no signal (fig. 2D and F). 
These results are consistent with a conserved function of piwi- 
and vasa-reldLted genes in germ cell maintenance and estab- 
lish their usefulness as germ cell markers in the leech. 

Temporal and Spatial Patterns of Hau-piwi and 
Hau-msa Expression in Development 
Leech development from deposition of the zygote (which is 
fertilized internally and then arrests at metaphase of meiosis I 
until deposition) to the juvenile has been divided into 1 1 stages 
(Fernandez and Stent 1982). To obtain an overview of Hau- 
piwil and Hau-vasal expression, we performed reverse tran- 
scription polymerase chain reaction (RT-PCR) for each devel- 
opmental stage as well as for the oocytes (fig. 3A and B). Both 
transcripts were detected at low levels in oocytes dissected 
from gravid individuals (stage 0) and in the final stages of 
development before the feeding juvenile (stages 9-11). 
Between these points, we observed two peaks of expression 
(fig. 3A and B). The first peak was associated with the zygote for 
Hau-piwil and with early cleavage in the case of Hau-vasal. 
Both transcripts declined to low levels during third cleavage 
(stage 4b) and then peaked again during stages 7-8. This 
second peak correlated with the production of segmental 
founder cells (blast cells) in sequential, anteroposterior pro- 
gression by the five bilateral pairs of lineage-restricted segmen- 
tation stem cells (teloblasts) that constitute a posterior growth 
zone in the embryos of clitellate annelids as described later. 

We have characterized the temporal and spatial expression 
of both duplicates of the piwi and vasa genes by ISH and 
found no significant differences between the expression 
of Hau-piwil and Hau-piwi2, or between Hau-vasal and 
Hau-vasa2. For brevity, only the results for Hau-piwil and 
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Fig. 1. Molecular phylogenies for piwi and wasa genes. Maximum likelihood trees for (A) piwi and (B) vasa family genes in relation to the respective 
outgroups plW and argonaut Helobdella robusta genes are indicated by shading. Numbers indicate bootstrap values equal to or greater than 50%. See 
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Fig. 2. HaU'piwil and Hau-vasal are expressed in gonads of the adult. Horizontal cryosections showing midbody segments of young adult Helobdella, 
stained with toluidine blue (A, B) or processed by ISH for Hau-piwil (C, D) or Hau-vasal (£, F). Anterior is up in all panels. For panels A, C, and £, the right 
side of each section is replaced by a line drawing indicating the reproductive (ovary, testisac) tissues. (A) Helobdella has one pair of ovaries, associated 
with segment M6, and six pairs of testisacs (five of which are visible here), associated with segments M8-M13. (B) Higher magnification views showing 
the arrangement of cells within a typical ovary and testisac. (C-F) Equivalent views of specimens processed by ISH show that Hau-piwil and Hau-vasal 
are expressed in both tissues. Boxes in (C) and (£) indicate areas shown at higher magnification in (D) and (F), respectively. Scale bar, 200 \xm in (A), (C), 
(£); 50|im in (B), (D), and (F). 



HaU'Vasal are presented below, but the following description 
applies to both piwi and both \/asa paralogs. 

Hau-piwil and Hau-vasal exhibited similar expression pat- 
terns at all developmental stages examined, from zygote to 
juvenile (stages 1-11). In the early zygote, both transcripts 
were confined to the region of the female pronucleus at the 
animal pole (fig. 3C). In the late zygote, the transcripts were 
associated with yolk-deficient domains of cytoplasm (telo- 
plasm) at the animal and vegetal poles of the embryo 
(fig. 3C). Teloplasm is enriched in mRNAs and mitochondria 
segregated to the D quadrant by unequal cleavage (Fernandez 
and Olea 1982; Fernandez et al. 1987; Holton et al. 1989, 1994). 
It is required for D quadrant determination (Astrow et al. 
1987). Hau-piwil and Hau-vasal mRNAs remained associated 
with teloplasm during cleavage but also appeared within the 
perinuclear cytoplasm in the macromeres. 

Stage 4b begins with the division of macromere (eight 
cell stage) to form blastomeres DM and DNOPQ, the precur- 
sors of mesodermal and ectodermal teloblasts, respectively. 
Both DM and DNOPQ inherit teloplasm. As in stage 3, both 
Hau-piwil and Hau-vasal transcripts at this stage were pre- 
sent in association with the teloplasm (of DM and DNOPQ) 
and the nucleoplasm of macromeres (fig. 3C). 

During blast cell (stages 6 and 7) and germinal band 
(stage 8) formation, transcripts were observed in all five telo- 
blast lineages (mesodermal and ectodermal). The expression 
of both transcripts was strong and homogeneous throughout 
the bandlets during these stages (figs. 3C and 4). During 
stage 9, the levels of both vasa and piwi transcripts decreased 
in anteroposterior progression (fig. 4). The residual expression 
was not uniform and appeared as transverse stripes on either 



side of the midline in each segment. During stage 10, these 
stripes of expression largely disappeared as well, leaving a 
single bilateral pair of Hau-piwil and Hau-\/asa1 positive 
spots in midbody segment 6 (segment M6), which we take 
to be the PGCs of the female germ line (figs. 4 and 5). Double- 
fluorescent ISH experiments confirmed that Hau-piwi and 
Hau-vasa were coexpressed in these cells (fig. 5). Staining in 
these two clusters of cells persisted through late stage 11, by 
which time the arms of the ovaries were visible in cross sec- 
tion (supplementary fig. SI, Supplementary Material online). 

No expression of piwi or vasa genes was observed in the 
male primordia of segments M8-M18 at any of the develop- 
mental stages examined. Curiously, we have previously ob- 
served robust expression of another germline marker, namely 
a nanos-class gene, in presumptive male germ line primordia 
during stages 9 and 10, well before the nanos homolog was 
detected in the presumptive female germ line (Kang et al. 
2002). From these two sets of observations, it would appear 
that in Helobdella, piwi and \/asa are expressed preferentially 
in the developing female germline, while nanos is expressed 
preferentially in the male germline. We have confirmed this 
conclusion by double ISH for Hau-nanos and Hau-piwi. 
Consistent with our previous results, Hau-nos is expressed 
preferentially in the presumptive male germline precursors 
in early stage 11, with the difference that we observed only 
ten pairs of spots (segments M8-M17) in H. austinensis com- 
pared with 11 (segments M8-M18) in H. robusta (Kang et al. 
2002; fig. 6A); by late stage 11, Hau-nos is coexpressed with 
Hau-piwi (fig. 6B) and Hau-vasa in the female germline pre- 
cursors and has disappeared from the presumptive male 
germline precursors (fig. 6B). 
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Fig. 3. Two peaks of Hau-piwi and Hau-wasa expression during early development. (A) Ethidium bromide-stained gel showing one of three replicate 
developmental RT-PCR experiments for Hau-piwil (30 cycles), Hau-vasal (30 cycles), and Hau-ISS-rRNA (25 cycles). (B) Graph showing the average 
amount (±SD) of Hau-piwil (filled circles) and Hau-vasa (open circles) relative to Hau-ISS-rRNA for the selected stages. (C) Representative embryos at 
selected stages that had been processed by ISH for Hro-piwil (top row) or Hro-wasal (bottom row): early stage 1 embryos (viewed obliquely) show 
transcript around the female pronucleus at the animal pole; stages late 1 and 2 (lateral views), 3 (animal pole view), and 4b (lateral view) show the 
association of both transcripts with teloplasm during cleavage (cells DNOPQ and DM are labeled); stage 7 (viewed obliqely with teloblasts and bandlets 
visible in dorsal posterior territory) show transcripts over nuclei (arrows), suggestive of active transcription; midstage 8 (viewed ventrally) shows broad 
accumulation of transcript throughout the germinal bands (small arrows) and germinal plate (large arrows). These data show an early peak of transcript 
accumulation in teloplasm, indicative of late maternal or very early zygotic transcription followed by decay to a minimum at the onset of fourth 
cleavage (stage 4b) and a rise in transcript levels during germinal band formation. Scale bar, 100 |im. 



Primorida Female Germline Cells Arise from 
Segmental Mesoderm in Helobdella 
in leeches and other clitellate annelids, segmental mesoderm 
arises in a lineage-dependent manner by the interdigitation of 
stereotyped clones of segmental mesodermal (sm) blast cells, 
which arise from the stem cell divisions of M teloblasts 
(Weisblat and Shankland 1985; Goto et al. 1999). To test 
the hypothesis that the p/iv/'and-i/^jsa'positive, female germ- 
line precursors in Helobdella arise from segmental mesoderm, 
as was found previously for the nanos-positwe presumptive 
male germline precursors (Kang et al. 2002), we carried out 
fluorescent ISH for piwi and/or \/oso on embryos in which M 
teloblasts had been injected with lineage tracers at various 
times. 

In one series of experiments, designed to test the hypoth- 
esis that the M teloblast lineage gives rise to the female germ- 
line, the newly formed left mesodermal teloblast (Ml) was 
injected with FarRed dextran amine at stage 5, and the resul- 
tant embryos were examined by confocal microscopy at early 
stage 1 1, after processing by ISH with probes for piwi and vasa: 
the left cluster of p/M//-and-i/asa-positive cells in segment M6 
was labeled with lineage tracer (supplementary fig. S2, 
Supplementary Material online), indicating that the female 
germline originates from the ipsilateral M lineage. 



In a second series of experiments, designed to identify 
which M blast cell gives rise to the female germline, the left 
mesodermal teloblast (MJ was injected with fluorescein dex- 
tran amine (FDA) 15 cell cycles after its birth; by this point, 
the M teloblasts had each generated six early mesodermal 
blast cells (em1-em6), which contribute to nonsegmental 
tissues (Gline et al. 2011) and then nine sm blast cells 
(sm1-sm9), whose clones are located sequentially in the 
four rostral segments (R1-R4, associated with the anterior 
sucker) and the first five midbody segments (M1-M5). The 
embryos were examined by confocal microscopy at stage 11, 
after processing by ISH with probes for piwi: the clusters of 
p/M//-expressing cells in segment M6 were labeled with lineage 
tracer on both sides in these embryos (fig. 7). On the left side, 
the p/u//-positive cluster was in the first labeled segment, in- 
dicating that the female PGGs originate from blast cell smIO 
in the M lineage. 

Discussion 

Differential Expression of Conserved Germ Line 
Markers in Developing Male and Female PGCs 
In the the work presented here, we have identified homologs 
of piwi and \/asa, evolutionarily conserved germ line markers, 
in the genome of the leech H. mbusta, and characterized their 
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Fig. 4. Restriction of Hau-piwil (left) and Hau-vasal (right) transcripts to female PGCs during later development. For each gene, the left: column shows 
side views of whole embryos at stages 9-11, anterior to left, while the right column shows ventral views of 5-10 segments, anterior up. The bottom 
image in each column shows a ventral view of an intact juvenile (stage latell, anterior to left). Black arrows indicate the prospective female PGCs 
(segment M6) in selected images; white arrows indicate the corresponding position in early stage 9 embryos, at which time differential staining of the 
PGCs is not obvious. Scale bar, 100 |im. 



expression in embryos of the closely related H. austinensis. The 
piwi genes from both species group in distinct piwil and piwil 
clades with pairs of piwi genes from other species, consistent 
with a duplication of the piwi gene early in bilaterian evolu- 
tion. By contrast, the two i^asa genes in Helobdella appear to 
represent a more recent duplication not seen in two other 
available lophotrochozoan genome sequences, from the mol- 
lusc L gigantea and the polychaete annelid C teleta. 

Leeches have discrete male and female reproductive or- 
gans, in contrast to the hybrid gonad (ovotestis) seen in most 
hermaphrodites; in Helobdella, a bilateral pair of ovaries is 
associated with midbody segment M6, and six pairs of testi- 
sacs lie in midbody segments M8-M13. In the adult, Hau-piwi 
and HaU'i^asa are coexpressed in the gonads and testisacs as 
expected (Mochizuki et al. 2001; Extavour and Akam 2003; 
Oyama and Shimizu 2007; Rebscher et al. 2007; Dill and Seaver 
2008; Kranz et al. 2010). Moreover, both piwi genes and both 
\/asa genes were expressed in both the male and female germ 
line in the adult. Thus, we found no evidence to suggest that 
the duplicated genes had evolved sex-specific functions. 

However, a striking result from our present work (summa- 
rized in fig. 8) is that, as levels of piwi and \/asa decline during 



stage 9, their transcripts become restricted specifically to the 
female germ cells (segment M6), whereas nanos transcripts 
become restricted first to presumptive male germ line pri- 
mordia during this same period (Kang et al. 2002). Thus, the 
differences do not arise simply from differences in the timing 
of the differentiation of the male and female organs but re- 
flect differential expression of nanos in developing male and 
\/asa and piwi in developing female germline. Intriguingly, 
these temporal differences in the expression of conserved 
germ line markers between the male and female germ cells 
correlate with greater developmental plasticity in male com- 
pared with the female germ line. First, as discussed previously, 
only six of the ten segments in which pairs of presumptive 
male germ line primordia express Hau-nos form definitive 
testisacs. Second, consistent with the predicted effects of 
social group size on sperm competition (Parker 1998), 
Helobdella exhibits postembryonic developmental adjust- 
ments in testisac volume, but not egg volume, depending 
on the number of animals in a mating group (Tan et al. 2004). 

We are unaware of observations regarding the differential 
expression of germline markers in other hermaphrodites that 
form separate ovaries and testes as does leech. In the ovotestis 
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Fig. 5. Coexpression of Hau-piwil and Hau-vasal in female PGCs. Two color fluorescent in situs for Hau-piwil (green; A, D, D') and Hau-vasal 
(red; B, E, £') and merged views (C, F, F'). Nuclei were counterstained with DAPI (blue). (A, B, C) Lateral view, anterior to left, showing the entire early 
stage 11 embryo. Labels indicate rostral segment 4 (R4), and midbody segments 1, 6, and 21 (Ml, M6, M21). (D, £, F) Ventral view, anterior up, showing 
midbody segments 5, 6, and 7 (M5, M6, M7). (D', E', F') Transverse section, ventral down, through midbody segment 6; labels in (D') indicate yolk 
(y) within the developing midgut (dotted contour) and the segmental ganglion (g). Scales bars, 100 |im in (A-C); 50|im in(D-F) and (D'-F'). 



of the nematode Caenorhabditis elegans, however, the piwi 
homolog prg'1 is required for spermatogenesis, and PRG-I is 
preferentially expressed in the male portion of the ovotestis 
during the spermatogenesis phase (Wang and Reinke 2008). 
C. elegans has also undergone a marked expansion within the 
argonaut gene family, which contains more than 20 argonaut- 
family genes. Two of these genes, alg'3 and alg'4, are required 
semiredundantly for spermatogenesis and ALG-B at least is 
preferentially expressed in the male portion of the ovotestis 
(Gonine et al. 2010). 

Delayed Segregation of Germ Line Precursors 
Gombined lineage tracing and situ hybridization shows that 
the presumptive female PGGs in Helobdella arise from the 
mesodermal lineage, specifically from the same M teloblast- 
derived blast cell that contributes most of the segmental 
mesoderm to midbody segment 6 (M6). We have previously 
shown that the presumptive male PGGs also arise from me- 
sodermal blast cells (Kang et al. 2002). Thus, as summarized in 



figure 8, the male and female germ line fates in Helobdella 
fates separate from one another only after 19 rounds of 
zygotic mitoses (and from somatic lineages even later), 
much later than in other invertebrate models. For example, 
germline is set aside at third cleavage in the amphipod crus- 
tacean Parhyale hawaiensis (Gerberding et al. 2002), at fourth 
cleavage in the nematode C elegans (Sulston et al. 1983), after 
nuclear division cycle eight in Dmsophila melanogaster (Foe 
and Alberts 1983), at the eighth and ninth divisions in the 
polychaete annelid P. dumerilii (Rebscher et al 2007), and at 
the 11th division in the mollusc Crepidula fornicata (Lyons 
et al. 2012). 

Leeches, together with earthworms and their relatives, 
comprise a monophyletic annelid subgroup called Glitellata 
(Erseus and Kallersjo 2004; Struck et al. 2011). The PGGs in 
another clitellate annelid, the oligochaete Tubifex, also arise 
from the M teloblast lineages (Goto et al. 1999), but whether 
from segmental blast cells remains to be determined. 
Whether or not inductive signals are required to specify the 
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Fig. 6. Differential expression of Hau-nanos and Hau-piwil in the male and female PGCs, respectively. Two color fluorescent in situs for Hau-piwil 
(green) and Hau-nanos (red). (A) Lateral views of early stage 11 embryo; Hflu-piwil is expressed in the female PGCs of segment M6 (top), Hau-nanos is 
expressed in male PGCs of segments M8-M17 (center) and there is no overlapping expression (bottom). (B) Ventral views of late stage 11 embryo; by 
this time, Hau-piwil (left) and Hau-nanos (center) are coexpressed in the female PGCs (right); Hau-nanos expression is not detected in the male 
germline at this stage. Scale bars, 100 |im in (A); 200 |im in (B). 



PGCs during embryogenesis in these species remains to be 
determined, but restoration of germline following amputa- 
tion in the oligochaete Enchytraeus (Sugio et al. 2008) must 
involve inductive signals of some sort. Together with the 
delayed segregation of germ line from somatic lineages in 
Helobdella, these results are consistent with the idea that 
some sort of epigenetic germ line specification was ancestral 
for this group (Extavour 2007). 

Broad Early Expression of Germ Line Markers 
The Helobdella piwi and \/asa genes are broadly expressed 
during early development (summarized in fig. 8) and then 
are gradually restricted to presumptive germ line cells (see 
next section). We have previously observed a similar pattern 
of expression for Hm-nanos (Kang et al. 2002), and similar 
patterns of broad early expression becoming gradually re- 
stricted to the PGCs have been observed for one or more 
of these conserved germ line markers for other annelid 
genera, including Tubifex (Oyama and Shimizu 2007), 
Platynereis (Rebscher et al. 2007), Capitella (Dill and Seaver 
2008; Giani et al. 2011), and Enchytraeus (Sugio et al. 2008). 
For Tubifex (like Helobdella, a clitellate annelid), it has been 
shown that restriction of i^asa to the germ line represents 
primarily differential stabilization of preexisting transcripts 
(Oyama and Shimizu 2007); by contrast, our RT-PCR results 



show that in H. austinensis both piwi and ifasa transcripts 
increase from the oocyte to the zygote or early cleavage 
stages, decline during mid-to-late cleavage, and then rise 
again during formation of the germinal bands and early ger- 
minal plate, indicative of a strong transcriptional component 
to their expression. 

The early expression of piwi, \jasa, and mms homologs 
throughout the embryo is seen not only in annelids but also 
in mollusks (Swartz et al. 2008; Kranz et al. 2010), cnidarians 
(Extavour et al. 2005), flatworms (Shibata et al. 1999; Pfister 
et al. 2008), arthropods (Nakao et al. 2006; Mito et al. 2008; 
Ozhan-Kizil et al. 2009), and echinoderms (Juliano et al. 2010). 
In several cases, a functional requirement for these genes in 
early development distinct from germ line specification has 
also been demonstrated. In H. robusta for example, a partial 
antisense morpholino knockdown of mms expression dis- 
rupts micromere positioning, germinal band elongation, and 
epiboly (Agee et al. 2006). It has been proposed that the broad 
early expression of piwi, msa, and mms homologs in both 
somatic and germline stem cells in flatworm and polychaete 
reflects a common origin for these cells during embryogenesis 
and that this may be the ancestral situation for metazoans 
(Rebscher et al. 2007). Elaboration of these ideas has led to the 
hypothesis of a broadly conserved germline multipotency 
program (GAAP) involving extensive translational repression 
(Juliano et al. 2010), which functions not only in PGCs but also 
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Fig. 7. Female PGCs originate from segmental mesodermal blast cell sm10. (A) Fluorescence micrographs of embryos in which the left: mesodermal 
teloblast (MJ was injected with FDA (green) 15 cell cycles after its birth, i.e., after the M teloblasts had each generated six presegmental blast cells (em1- 
em6) and nine segmental blast cells (sm1-sm9). The embryos were examined by confocal microscopy at stage 11, after processing by ISH with probes 
for HaU'piwH (red). The p/M//-expressing cell(s) in segment M6 were labeled with lineage tracer on the left side. (A) Lateral view showing left side, anterior 
to left. (B) Ventral view, anterior up showing roughly eight segments. Dotted contour outlines the visible portion of the embryo; on the left (apparent 
right) side, the p/n//- positive cluster is in the first labeled segment, indicating that the female PGCs originate from blast cell smIO in the M lineage. 
(C) Transverse section of a similar embryo. Coarse dotted contour outlines the embryo; finer contours outline the yolk (y) and the ventral ganglion (g), 
Scale bar, 100 |im in (A); 50|im in (B); 30|im in (C). 



in an evolutionarily ancient type of PriSCs (Solana 201 3), from 
which germ cells, somatic stem cells (neoblasts), and somatic 
tissues may arise. Our present results support these ideas. 

Speculation on a Possible Developmental Basis for 
Acceleration of Genome Rearrangements in 
Helobdella Relative to Other Animals 
Recent analyses of three lophotrochozoan genomes (Simakov 
et al. 2013) revealed a dramatic acceleration of gene rearran- 
gements relative to protein evolution in the leech H. robusta 
that was not seen in the other species sequenced. Our present 
results suggest an interplay of developmental and reproduc- 
tive strategies that might explain this. 

Hermaphroditism is an ancestral trait for clitellate annelids. 
It is commonly held that clitellates are incapable of self-fertil- 
ization, but self-fertility has been observed for H. triserialis 



(Wedeen et al. 1990), H. robusta (Weisblat DA, unpublished 
observations), and H. papillornata (Tan et al. 2004), repre- 
senting at least two distinct clades within the genus Helobdella 
(Siddall and Borda 2003; Bely and Weisblat 2006). Helobdella 
austinensis (Kutschera et al. 2013) is in the H. robusta clade 
(Bely and Weisblat 2006) but is not self-fertile (Levine M, 
Weisblat DA, unpublished data); given the apparent preva- 
lence of self-fertility in the Helobdella genus, the self-infertility 
in H. austinensis may represent a secondary loss. Self-fertiliza- 
tion, an extreme form of inbreeding, provides a partial expla- 
nation for the accelerated genome rearrangements: otherwise 
lethal genetic translocations that move genes from one chro- 
mosome (linkage group) to another can be rescued by com- 
bining gametes that have the same translocation. 

But this explanation alone is unsatisfactory. The nematode 
C elegans is also a self-fertile hermaphrodite and has a much 
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DM "I dm" 
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Stage 2 



Stage 4b 




(animal pole view) (animal pole view) 



Early stage 11 
(dorsal view) 



Late stage 11 
(dorsal view) 



Late stage 8 
(ventral view) 



Adult 
(dorsal view) 






• Hau-piwi and vasa # Hau-nanos • Hau-piwi, vasa, and nanos 

Fig. 8. Diagrammatic summary of germ line marker expression and PGC segregation in Helobdella. Lineage diagram depicting cell divisions leading from 
the zygote to the segmental founder cells from which female and male PGCs arise, with line drawings showing the expression of conserved germ line 
markers at selected stages and in the adult leech (nanos, blue; coexpression of piwi and vasa red; coexpression of nanos, piwi and vasa, purple). Time line 
is approximate: from zygote to stage 7 takes about 54 h at 23 °C; teloblast cell cycles are 1-2 h. Bilateral M teloblasts arise at seventh cleavage from the 
division of DM'' (micromere 4d in standard nomenclature for spiralian embryos); each M teloblast undergoes six further rounds of stem cell division 
that produce nonsegmental early mesoderm blast cells (em1-em6) before producing its first segmental blast cell (sm1) (Gline et al. 2011); finally, six 
further rounds of mitosis occur before the precursor of the female germ line (sm6) are separated from the precursor of the male germ line (which arise 
from the M teloblast starting with sm8). As shown in figure 3C, coexpression of all three markers is prominent in teloplasm and teloblast precursors (but 
is also seen in other cells) during cleavage stages (e.g., stages 2 and 4b), and in the germinal bands and early germinal plate (stage 8). By early stage 11, 
expression of piwi and vasa has become restricted to the females PGCs (segment M6), while nanos is restricted to prospective male PGCs (segments 
MS-MI 7 in H. austinensis; and also segment Ml 8 in K wbusta, indicated by paler shading). By late stage 11, expression of all three genes is restricted to 
the female PGCs. The point at which expression resumes in adult male PGCs remains to be determined. See text for further details. In the adult, ovaries 
are associated with segment M6 and definitive testisacs with segments M8-M13. 



faster generation time (4 days for C elegans vs. 9-10 weeks for 
H. wbusta under laboratory conditions). Accordingly, its 
genome has diverged much more rapidly from the inferred 
ancestral protostome than has H. wbusta by most criteria 
except that of genome rearrangements; traces of macrosyn- 
teny are still evident in the genome of C elegans but not in 
that of Helobdella (Simakov et al. 2013). 

Our analysis of germline formation in Helobdella suggests a 
plausible resolution to this paradox. Whereas Caenorhabditis 
segregates germline from somatic cells at the fourth zygotic 
mitosis (and generates male and female gametes by a some- 
what indeterminate combination of proliferative and stem 
cell divisions thereafter; Kimble and Hirsch 1979), ipsilateral 
male and female PGC lineages in Helobdella are not segre- 
gated from each other until 19 rounds of zygotic mitosis have 
occurred (fig. 8). A priori, this extensively shared lineage in- 
creases the opportunity for genetic translocations to occur in 
cells that will contribute both male and female PGCs and 



which therefore might be rescued by self-fertilization. 
Conversely, self-fertilization selects against any homozygous 
lethal mutations that arise in a common precursor of male 
and female PGCs. Thus, the combination of reproductive 
strategy (self-fertile hermaphroditism) and developmental 
features (delayed segregation of male and female PGCs) 
seen in modern H. wbusta suggests an explanation for how 
the lineage from which this species arose might have accu- 
mulated extensive gene rearrangements, while other aspects 
of its genome relatively remained stable. 

Testing this explanation for the extensive rearrangements 
of the H. wbusta genome relative to other animals is not 
possible at present, given the scant information regarding 
genome organization and reproductive strategies among re- 
lated taxa. Similar challenges arise in interpreting results from 
other examples of self-fertilizing hermaphrodites. The killifish 
Kryptolebias marmomtus (Teleostei: Cyprinodontiformes: 
Rivulidae) and a sibling species K. ocellatus are the only 
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vertebrate hermaphrodites for which self-fertilization has 
been well documented (Tatarenkov et al. 2009). Microsatellite 
analysis of four Kryptolebias species indicates that the capacity 
for self-fertilization in K. marmoratus and K. ocellatus has been 
maintained for several hundred thousand years or more 
(Tatarenkov et al. 2009); analysis of locally most common 
multilocus genotypes suggests that fertilization assurance 
rather than mulitlocus coadaptation is the driving force fa- 
voring self-fertility in this species (Avise and Tatarenkov 2012). 
Kryptolebias marmoratus shows chromosomal constancy 
over a large geographic range, but whether this represents 
recent, rapid spread or selection against the expected rear- 
rangements cannot be determined (Sola et al. 1997). 

In the present case, support for our explanation would be 
provided by finding that Helobdella is the only clitellate group 
to have evolved self-fertilization and is the only species in this 
group with such extensive genome rearrangements. Alterna- 
tively, if genomic rearrangments are a characteristic of clitel- 
late annelids as a group, it could reflect earlier gains and losses 
of self-fertilization within the clitellate lineage, or simply more 
extensive inbreeding associated with fluctuating populations 
of hermaphroditic animals in ephemeral freshwater habitats. 

Materials and Methods 

Animals 

A laboratory colony of the species newly described as 
H. austinensis (Kutschera et al. 2013) was raised and embryos 
were collected as described elsewhere (Gline et al. 2011). 

Gene Identification, Cloning, and In Vitro 
Transcription 

Putative PAZ and piwi domains (for Piwi and Argonaut) and 
DEAD/DEAH box helicase domains (for vasa and PL10) 
encoded by the annotated H. robusta genome (DOE, Joint 
Genome Institute) were identified by a Pfam search, using 
accession numbers PF02170, PF02171, and PF00270, respec- 
tively. As a further step, gene models of piwi/Ago and vasa/ 
PL10 orthologs were used to perform BLASTP similarity 
search of the H. robusta genome. For each locus, the sets of 
gene models were assembled into a single predicted best-fit 
model for each gene using the SeqManll software (DNASTAR, 
Madison, Wl). 

Specific primers (piwil forward: ACCTCCATACTGGAGTA 
TGATGAGGGAATA; piwil reverse: TAGTTTGTTGGCCAGAG 
TGTTACTAGGTGA; \/asa1 forward: TTCCTGCTTCCTGTGTT 
GACGGGTATGTTG; \/asa1 reverse: GCCACATCTACCAGTTC 
TACCAATCCTATG) were designed against the piwil and 
i^asal gene models and amplified from cDNA prepared 
from the closely related H. austinensis. PGR products were 
gel extracted with a Qiagen Gel Extraction Kit and ligated 
into pGEM-T easy vector according to manusfacturer's in- 
structions. These clones were introduced into Escherichia 
coli DH105 competent cells by heat shock, plated onto 
solid media (15 g agarose per liter LB broth plus 50|ig/ml 
ampicillin) and grown overnight at 37 °G. Golonies were se- 
lected and grown overnight in LB broth with 1% ampicillin. 
Plasm ids were isolated using the Qiagen Miniprep Kit 



following manufacturer's instructions and sequenced com- 
mercially on an ABI platform. The MEGAscript (Ambion) 
kit was used to make digoxygenin (Dig-11-UTP) and/or 
biotin (Biotin-16-UTP) riboprobe according to manufac- 
turer's instructions. 

Alignment and Phylogenetic Analyses 

Piwi/Argonaut 

A data set was built with orthologs of piwi and argonaut 
from representatives of the major metazoan clades with a 
total of 28 amino acid sequences (supplementary tables SI 
and S2, Supplementary Material online). These sequences 
were prealigned using AMAP (Larkin et al. 2007; Schwartz 
and Pachter 2007) followed by manual editing in Se-AI 
v2.0a11 (http://tree.bio.ed.ac.uk/software/seal/, last accessed 
October 28, 2013). Gblocks were used to exclude poorly 
aligned amino acid positions under the least stringent param- 
eters as implemented in the online version of the program 
(Gastresana 2000). The resulting alignment comprised 512 
amino acid positions (TREEBASE 14911 [http://purl.org/ 
phylo/treebase/phylows/study/TB2:S1491 1, last accessed 
November 8, 2013]). The Whelan and Goldman (WAG) 
model of protein evolution (Whelan and Goldman 2001) 
with gamma distribution + invariant sites was implemented 
in all analyses as selected by ProtTest v1.4 as the model that 
best fits the data (Abascal et al. 2005). Maximum likelihood 
analyses were executed in RAxML-v7.0.4 (Stamatakis et al. 
2005), for which support was obtained from 1,000 bootstrap 
replicates. Trees were viewed and edited in Dendroscope 
V2.2.2 (Huson et al. 2007). 

Vasa/PUO 

A matrix comprising 26 Vasa and PL10 orthologs representing 
the major metazoan clades was constructed. This data set 
(comprising 370 amino acid residues) was analyzed as for 
Piwi/Argonaut (see earlier). 

Developmental RT-PCR 

Total RNA samples were prepared from H. austinensis em- 
bryos at selected stages with RNA Wiz (Ambion) according to 
the manufacturer's instructions, using 40 embryos for each 
stage. The total RNA obtained (3 |ig) was reverse transcribed 
using a first-strand cDNA synthesis kit (BD Biosciences, Palo 
Alto, GA). The PGR mixture (50 |il) contained lOx Taq buffer, 
0.3 U Taq polymerase (Perkin-Elmer), 2.5 |iM of dNTPs, 
5 pmol of each set of primers, and 50 ng of cDNA from the 
selected stages as template (fig. 2A). To amplify Hau-piwil- 
and /-/c?t/-i/asc? 7 -specific fragments, a pair of PGR primers was 
designed as follows: Piwil-F: 5^-GGGGATATTGGGAAGATAAG 
GGTGAGAAA-3^ Piwil -R: 5'-TATTGGGTGATGATAGTGGAG 
TATGGAGGT-3^ Vasal -F: 5^-TTGGTGGTTGGTGTGTT GAG 
GGGTATGTTG-3^ Vasal R: 5^-GGGAGATGTAGGAGTTGTAG 
GATTGGTATG-3^ The PGR reactions were performed under 
the following cycling conditions: an initial denaturation at 
94 °G for 5 min followed by 25-35 cycles of denaturation at 
94 °G for 30 s, annealing at 60-65 °G (depending on the 
primer set) for 30 s, and elongation at 72 °G for 1 min, with 
a final elongation step at 72 °G for 10 min. A 10-|il aliquot of 
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each PCR reaction was removed after 25 cycles, while the 
remaining material underwent five more cycles of amplifica- 
tion. Each primer set yielded a PCR product of 387 and 882 bp 
in length for Hau-piwil and HaU'\/asa1, respectively. The 18S 
rRNA sequence was used as an internal standard 
(QuantumRNA, Ambion). The extent of amplification was 
chosen empirically to avoid saturation of the amplified 
bands. To quantify PCR products, each sample was run in a 
1.5% agarose gel and stained with ethidium bromide. Band 
intensity was measured with an Alphaimager (Alpha 
Innotech Corp.) using the Alphaease (v3.3b) program. 

Histology 

To visualize differentiated testisacs and ovaries, young adult 
leeches were fixed in 4% formaldehyde in PBS overnight at 
4°C and mounted, sectioned, and dried as described by Cho 
et al. (2010). Cryosections were mounted on glass slides and 
stained with freshly filtered toluidine blue (1 mg/ml 2.5% 
aqueous sodium carbonate, pH 11.1) before placing the 
coverslip. 

ISH Lineage Tracing and Cross Sectioning 
All ISHs, double-fluorescent ISHs, lineage tracing, and cross 
sectioning were performed as previously described (Cho et al. 
2010). Embryos were stained with 4^6-diamidino-2-phenylin- 
dole (DAPI, Sigma) to visualize cell nuclei. Probe lengths were 
1,848 bp for Hau-piwil, and 876 bp for HaU'i^asal. 

ISH on Cryosections for Adult Leech 
Adult leeches were fixed at 4°C overnight (4% paraformal- 
dehyde in 1 X PBS) and dehydrated in methanol series, equil- 
ibrated in 30% sucrose and mounted in OCT (Tissue-Tek) and 
stored frozen before sectioning on a cryostat (8-|im sections). 
Prepared sections were dried (50 °C, 30min) and stored at 
— 70°C until further use. The specimens were incubated at 
37 °C for 10 min in 20 |ig/ml pronase E in PBT (PBS plus 0.1% 
Tween20, pH 7.4), then rinsed three times in PBT and incu- 
bated in 0.3% H2O2 in methanol at room temperature 
for 20 min. Sections were then postfixed for 20 min in 4% 
paraformaldehyde and finally rinsed twice in PBT. 
Prehybridization was carried out at 67 °C for 30 min in hy- 
bridization buffer (50% deionized formamide, 5x SSC, 1x 
Denhardt's solution, 0.1% CHAPS, 100|ig/ml heparin, 0.1% 
Tween 20, 100|ig/ml tRNA). The prehybridization buffer 
was replaced with fresh hybridization buffer containing 
2ng/ml of the corresponding probe (Hau-piwH or Hau- 
\/asa1) and specimens were incubated at 67 °C overnight. 
Washed specimens were incubated at room temperature 
for 1.5 h in 1% blocking reagents (Roche) in PBT, then incu- 
bated at 4°C for 16 h with 1/1,000 Anti-DIG/AP antibody 
(Roche) in 1% blocking reagents. After incubation, the color 
reaction was carried out using nitro blue tetrazolium chlo- 
ride/5-bromo-4-chloro-3-indoyl-phosphate (Roche) by stan- 
dard procedures. Stained specimens were dehydrated in 
ethanol, mounted in Fluoromount-G (Southern Biotech), 
and examined by brightfield microscopy. 



Supplementat7 Material 

Supplementary figures SI and S2 and table SI are available at 
Molecular Biology and Evolution online (http://www.mbe. 
oxfordjournals.org/). 
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